Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



The shape of the blue/UV continuum of B3-VLA radio 
quasars: Dependence on redshift, blue/UV luminosity and 
radio power 

Ruth Carballo 1 ' 2 , J.Ignacio Gonzalez-Serrano, 1 , Chris R. Benn 3 , 
Sebastian F. Sanchez 1 ' 2 and Mario Vigotti 4 

1 Instituto de Fisica de Cantabria (CSIC-Universidad de Cantabria). Facultad de Ciencias, 39005 Santander, Spain 

2 Departamento de Fisica Moderna, Universidad de Cantabria. Facultad de Ciencias, 39005 Santander, Spain 

3 Isaac Newton Group of Telescopes, Royal Greenwich Observatory, Apdo. de correos 321, 38780 Santa Cruz de La Palma, Tenerife, Spain 

4 Istituto di Radioastronomia, CNR, via Gobetti 101, 40129 Bologna, Italy 



1 February 2008 



ABSTRACT 

We present UBVR photometry of a sample of 73 radio quasars, about 80 per cent 
complete, with redshifts 0.4-2.8. From these data the shape of the spectral energy- 
distribution (SED) in the rest-frame blue/UV is analysed, using the individual sources 
as well as through broad-band composite SEDs. The SEDs of the individual sources are 
generally well fitted with power-laws, with slopes a ranging from 0.4 to -1.7 {S v oc v a ). 
Two sources with a < —1.6 were excluded from the general study for having very red 
SEDs, significantly deviating with respect to the remaining sources. The composite 
SEDs cover the range ~1300-4500 A and the only emission feature apparent from the 
broad-band spectra is the CIVA1549 line, in agreement with expectations from line 
equivalent width measuremets of radio-loud quasars from the literature. The shape of 
the composites in the logS^-log^ plane exhibits a break at around 3000 A where the 
spectrum changes from abiuc = 0.11 ± 0.16 at A > 3000 A to auv = —0.66 ± 0.15 at 
A < 3000 A. Although the broad-band spectral points are expected to include some 
masked contamination from emission lines/bumps, the break cannot be explained by 
line/bump emission, and most likely reflects an intrinsic trend in the continuum. 

The continuum shape is shown to depend on redshift. For the quasars with z < 1.2 
we find abiuc = 0.21±0.16 and arjv = — 0.87±0.20, i.e. a higher steepening. For z > 1.2 
auv is more flat, —0.48 ± 0.12, and there are too few spectral points longward of 3000 
A to obtain abiuc and analyse the presence of the 3000 A break. A trend similar to 
that between qvjv and z is found between auv and luminosity at 2400 A, £2400 j with 
luminous quasars exhibiting a harder spectrum. 

The data show an intrinsic correlation between L2400 and the radio power at 408 
MHz, not related to selection effects or independent cosmic evolution. The correlations 
Oi\jy-z, auv~^24oo and L^q-z appear to be consistent with accretion disc models 
with approximately constant black hole mass and accretion rates decreasing with time. 
If the trends L2A0Q-Z and P^-z are predominantly related to a selection bias, rather 
than cosmic evolution, only one of the correlations auv~^2400 or a uv~z need to be 
intrinsic. 
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1 INTRODUCTION 

The study of the shape of the optical/UV continuum of 
quasars is an essential tool to test the different emission 
models invoked to explain the observed radiation, and un- 
derstand their characteristic parameters. Whereas the over- 



all quasar spectrum from the infrared to the X-ray has ap- 
proximately a power-law form, S v oc v a , with spectral in- 
dex about — 1 , in the optical/UV region there is a bump 
on top of the power-law continuum known as the big blue 
bump (Elvis et al. 1987 & Sanders et al. 1989). The big blue 
bump is generally interpreted with a two component model, 
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consisting of the underlying power-law and black-body-like 
emission from an accretion disc (hereafter AD) around a 
massive black hole (Malkan 1983, Czerny & Elvis 1987). 
The determination of the shape and strength of the ioniz- 
ing UV continumm is essential to constrain the accretion 
disc parameters, for the modelling of the broad line region 
of quasars, and other aspects such as the determination of 
the fc-corrections for quasars, which affect the derivation of 
their luminosity functions. 

Empirically, the shape of the optical/UV continuum of 
quasars is generally parameterized by a power law, although 
this is a local approximation, with the overall shape being 
more complicated. Measurements of the optical/UV contin- 
uum shape at fixed rest-frame wavelengths have been ob- 
tained by O'Brien et al. (1988), and more recently by Na- 
tali et al. (1998), on the basis of low-resolution spectroscopic 
data. O'Brien et al. selected their quasars from availability of 
IUE spectra and published X-ray fluxes, and found a mean 
spectral index of -0.67 ± 0.05 for the range 1900-1215 A. 
The spectral index distribution of radio-loud and radio-quiet 
quasars appeared to be similar. The authors found a small 
hardening with redshift, with a ranging from —0.87 ± 0.07 
for z < 1.1 to —0.46 ± 0.05 for z > 1.1, and a trend with 
luminosity at 1450 A, in the sense that the more luminous 
quasars had harder spectra. From a joint regression analysis 
including spectral index, redshift and luminosity, the au- 
thors concluded that the dominant correlation was between 
q and z, and that the trend between a and luminosity was 
due to the correlation of both variables with redshift. Natali 
ct al. (1998) used a complete sample of optically selected 
radio quiet quasars, and found that the spectra in the range 
5500-1200 A showed an abrupt change around the so called 
"3000 A emission feaure" (Wills, Netzer & Wills 1985), with 
ol — 0.15 for A > 3000 A and a ~ -0.65 for A < 3000 A. 

Francis et al. (1991) analysed the composite spectrum 
obtained from optical spectra of around 700 quasars from the 
LBQS (Hewett et al. 1991). The individual spectra yielded 
slopes ranging from 1 to -1.5 and a median value a=-0.3 in 
the range 5800-1300 A. Francis et al. fitted to the composite 
spectrum a curved underlying continuum (a cubic spline), 
corresponding to a ~ 0.0 for A > 3000 A and a ~ -0.6 
below this limit. The authors noted that the high- luminosity 
quasars had harder spectra than the low luminosity ones. 

The samples studied by Francis et al. (1991) and Na- 
tali et al. (1998) select quasars on the basis of the presence 
of UV excess (although the LBQS includes additional in- 
dependent selection criteria), and are then biased against 
quasars with red colours. A similar bias is likely present in 
the data of O'Brien et al. (1988), who selected the quasars 
by availability of UV and X-ray measurements. In this pa- 
per we present a study of the shape of the blue/UV spec- 
trum of quasars from the B3-VLA Quasar Sample (Vigotti 
et al. 1997) selected in radio, at 408 MHz. Optical selection 
bisases are significantly lessened in this sample, allowing for 
a thorough investigation of the shape of their blue/UV spec- 
trum. Since the B3-VLA quasars have been selected at a low 
frequency, where steep-spectrum, extended emission domi- 
nates, the sample minimizes the inclusion of core-dominated 
quasars, for which an additional component of relativisti- 
cally beamed optical synchrotron is likely present (Browne 
& Wright 1985). 

The analysis of the blue/UV continuum of the B3- 



VLA quasars is based on UBVR photometry of around 
70 sources, with redshifts in the range 0.4-2.8. The quasar 
SEDs are studied both individually and from created com- 
posite spectra, and the dependence of the blue/UV slope 
and luminosity on redshift and radio properties is analysed. 

Near infrared K-banA imaging of 52 quasars in this work 
was presented by Carballo et al. (1998). For sixteen sources 
the images revealed extended emission, most likely related 
to starlight emission from the host galaxy. 



2 THE SAMPLE 

The B3 survey (Ficarra, Grueff & Tomasetti 1985) cata- 
logues sources to a radio flux-density limit 5*408 MHz = 0.1 
Jy. From the B3 Vigotti et al. (1989) selected 1050 radio 
sources, the B3-VLA Sample, consisting of five complete 
subsamples in the flux density ranges 0.1 — 0.2 Jy, 0.2 — 0.4 
Jy, 0.4 — 0.8 Jy, 0.8— 1.6 Jy and S408 >1.6 Jy, all mapped at 
the VLA. Candidate quasar identifications (objects of any 
colour appearing starlike to the eye) were sought on the 
POSS-I red plates down to the plate limit, R ~ 20, yield- 
ing a sample of 172 quasar candidates. Optical spectroscopy 
was obtained for all the candidates and 125 were confirmed 
as quasars, forming the B3-VLA Quasar Sample. The sam- 
ple covers the redshift range 2 = 0.4 — 2.8, with mean red- 
shift z = 1.16, and radio powers P408 MHz ~10 33 -10 36 erg 
s _1 Hz -1 (adopting # =50 km s^Mpc -1 and f2 = l). The 
sample of quasar candidates and the final B3-VLA Quasar 
Sample are described in Vigotti et al. (1997). 

The optical incompleteness of the radio quasar sam- 
ple, i.e. the fraction of quasars fainter than the optical limit 
of R ~ 20 mag, depends on the radio flux. From _R-band 
photometry of the complete radio quasar sample presented 
by Willot et al. (1998), we infer that the fraction of ra- 
dio quasars with R > 20 is 35 per cent for the flux range 
0.4 < S408 < 1-0 Jy (mean flux density 0.66 Jy) and 45 
per cent for 0.2 — 0.4 Jy (mean flux density 0.28 Jy). The 
average B — R colour of these quasars is 0.6 (similar to the 
value we found for the B3-VLA quasars on Section 4.2). For 
the quasars with S408 > 1 Jy we can use the distribution 
of S-band magnitudes obtained by Serjeant et al. (1998) 
for the Molonglo-APM Quasar Survey (average flux den- 
sity 5408 ~ 1.7 Jy). Adopting as limit B — 20.6, equivalent 
to R — 20 for typical radio quasar colours, the fraction of 
quasars with 54os > 1 Jy fainter than this limit is around 
15 per cent. 

The B3-VLA Quasar Sample contains 64 quasars with 
5 > 0.8 Jy and average flux density of 2.25 Jy, 32 with 
5 = 0.4 — 0.8 Jy and average flux density of 0.53 Jy and 29 
quasars with 5 = 0.1 — 0.4 Jy and average flux density of 0.22 
Jy. Adopting for the three groups optical incompleteness of 
15 per cent, 25 per cent and 45 per cent respectively, the 
estimated optical completeness for the total sample would 
be around 75 per cent, improving to 80 per cent for 5 > 0.4 

Jy- 

The present work is based on UBVR photometry of a 
representative group of 73 quasars from the B3-VLA Quasar 
Sample. The quasars were selected to have right ascensions 
in the range 7 -15 ft and comprise the 44 with S408 > 0.8 
Jy, the 23 with 0.4 < 5 4 os < 0.8 Jy, the 4 with 0.3 < 
5 40 8 < 0.4 and 2 out of 15 with 0.K 5 40 8 < 0.3. The 



© 0000 RAS, MNRAS 000, 000-000 



Blue/ UV continuum of B3-VLA radio quasars 3 



sample is equivalent to the Quasar Sample (except for the 
R.A. constraint) for S > 0.4 Jy, and includes only a few 
quasars with fainter radio fluxes, although generally close 
to the limit. We estimate therefore the optical completeness 
of the studied sample to be about 80 per cent. 



3 OBSERVATIONS, DATA REDUCTION AND 
PHOTOMETRIC CALIBRATION 

The optical images were obtained on 1997 February 5—8, 
using a 1024x1024 TEK CCD at the Cassegrain focus of 
the 1.0-m JKT on La Palma (Spain), and on March 11 — 12, 
using a 2048 x2048 SITe CCD at the Cassegrain focus of 
the 2.2-m telescope on Calar Alto (Spain). U,B,V and R 
standard Johnson filters were used and the pixel scale was 
0.33 arcsec pixel -1 for the JKT and 0.533 ascsec pixel -1 for 
the 2.2-m telescope. The field of view was ~6'x6' for all the 
images. A standard observing procedure was used. Bias and 
sky flat fields for each night and filter were obtained dur- 
ing the twilight. Faint photometric standards (Landolt 1992 
and references therein) were observed each night in order to 
obtain the flux calibration. The exposure time was different 
for each quasar and filter, ranging from 60 s to 1200 s, and 
was set according to the red and blue APM magnitudes from 
POSS-I. The seeing varied from ~ 1.8 arcsec, on February 
5, 6, 7 and March 11, to ~ 2.2 arcsec on February 8 and 
March 12. The data were reduced using standard tasks in 
the iraf software package Flat field correction was better 
than 0.5 per cent. Large exposure time images were cleaned 
from cosmic-ray hits automatically. 

Instrumental magnitudes for the standard stars were 
measured in circular apertures of 13 arcsec diameter. The 
flux calibration was obtained, as a first step, assuming no 
colour effects in any of the bands. For all the nights at both 
telescopes this assumption was proved to be correct for fil- 
ters B, V, and R. At the U band, for two of the nights, we 
took colour effects into account, introducing the colour term 
—k"(U — B). Table 1 lists the results of the photometric cal- 
ibration showing the extinction and colour coefficients for 
each night /filter together with the rms of the fits. 

All nights were photometric except the second half of 
February 5th, during which we had clouds. In this part of the 
night we only obtained data for 3 objects in two bands. The 
first part of the night was photometric and the calibration 
data listed in Table 1 correspond to this part of the night. 



4 OPTICAL PHOTOMETRY OF THE 
QUASARS 

4.1 UBVR magnitudes 

Quasar magnitudes were measured on the images using the 
same apertures as for the photometric standards. In one 
case (B3 0724+396) a nearby star was included within the 
aperture and it was subtracted by modelling it with a two- 
dimensional PSF. Measurement errors ranged from less than 
0.01 mag to 0.9 mag, and the typical values were lower than 

* IRAF is distributed by the NOAO, which is operated by AURA, 
Inc., under contract to the NSF 



0.15 mag (85 per cent of the data). Apparent UBVR mag- 
nitudes, not corrected from Galactic extinction, and their 
errors are listed in Table 2, along with the observing dates 
and the colour excess E(B — V) towards each object, ob- 
tained from the JV(Hi) maps by Burnstein & Heiles (1982). 
The redshift and radio power of the sources is also listed in 
the Table. 

Some of the quasars (11) were observed with the same 
filter twice on the same night or on the two runs. In these 
cases the difference in magnitudes never excedeed 0.2 mag, 
and the quoted magnitudes correspond to the average value. 
For three quasars (B3 0836+426, B3 0859+470 and B3 
0906+430) the U and V-band data were obtained the second 
part of Feb 5th, which was cloudy. We have made a crude 
estimate of these magnitudes on the basis of the photomet- 
ric standards measured over the same period, and assigned 
them an error of ~ 0.3 mag. These objects, however, have 
not been used for any further analysis. 

For the 70 quasars observed under good photomet- 
ric conditions rest-frame SEDs were built, using the zero- 
magnitude flux densities from Johnson (1966; U,B and R 
bands) and Wamsteker (1981; V band). For the Galactic ex- 
tinction correction Rieke & Lebofsky (1985) reddening law 
was used. The corrections were lower than 0.05 mag for 80 
per cent of the sources. For the remaining sources the cor- 
rections for the U band, which is the most affected, ranged 
from 0.05 to 0.56 mag, with a median value of 0.19. Figure 
1 shows the SEDs of the 70 quasars, plotted in order of in- 
creasing redshift. The SEDs will be always referred to the 
\ogS v -\ogv plane. 

Three sources in Table 2 (B3 0955+387, B3 1312+393 
and B3 1444+417, see also Fig. 1) have magnitude errors 
larger than 0.3 mags in several bands. In addition, three 
sources have abrupt changes in their SEDs (B3 1317+380, 
B3 0726+431 and B3 0922+425), probably related to intrin- 
sic variability or not-understood errors. These six sources 
will not be considered for the analysis presented in the forth- 
coming sections. 

Histograms of the U,B,V and R magnitudes (corrected 
for Galactic extinction) for the 64 quasars with good pho- 
tometry are presented on Figure 2. The ii-band histogram 
shows that most of the quasars have R magnitudes well be- 
low the POSS-I limit of R ~ 20 used for the quasar identifi- 
cations, confirming that this limit guaranties a rather high 
optical completeness. 

We have indicated on Fig. 1 the wavelengths of the 
strongest quasar emission lines in the studied range, like H/3, 
MgllA2798, Cm]A1909, CivA1549 and Lya. For each spec- 
tral point on the figure we have plotted the covered FWHM, 
assuming the standard FWHMs (observer frame) for the 
U,B,V and R bands from Johnson (1966). We infer from 
the figure that about 45 per cent of the spectral points do 
not include within half the filter bandpass the central wave- 
lengths of any of the emission lines listed above. Broader 
emission features, such as the Fell and Fell+Balmer emis- 
sion bumps in the ranges 2250-2650 A and 3100-3800 A, 
could also affect the measured broad-band fluxes. 

Table 3 presents average equivalent width (EW) mea- 
surements of these emission lines for several quasar samples 
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Table 1. Summary of photometric calibration 



Date 


k(U) 


k"(U - B) 


rms 


KB) 


rms 


k(V) 


rms 


k(R) 


rms 


5 Feb 


1.55 ±0.15 




0.024 






0.54 ±0.12 


0.009 






6 Feb 


0.34 ±0.07 




0.014 






0.12 ± 0.10 


0.027 






7 Feb 


0.54 ±0.10 




0.051 


0.24 ±0.02 


0.012 


0.21 ± 0.02 


0.020 


0.15 ±0.02 


0.015 


8 Feb 








0.28 ±0.04 


0.021 






0.11 ±0.01 


0.008 


11 Mar 


0.48 ± 0.05 


-0.250 ±0.023 


0.042 


0.19 ±0.03 


0.028 


0.15 ± 0.01 


0.013 


0.09 ±0.01 


0.008 


12 Mar 


0.34 ±0.05 


-0.163 ±0.030 


0.052 


0.14 ±0.03 


0.042 


0.13 ± 0.02 


0.027 


0.06 ± 0.02 


0.022 



Table 2. Photometry of the B3-VLA quasars 



B3 name 




ergs — ^ Hz — ^ 


u 


a(U) 


B 




v 




R 




Da/tc of obs. 
U/B/V/R 


E(B — V) 


0701±392 


1.283 


34.99 


18.50 


0.14 


20.11 


0.17 


19.69 


0.12 


18.89 


0.13 


6/8/6/8 


0.118 


0704+384 


0.579 


34.72 


18.82 


0.12 


19.08 


0.18 


18.84 


0.14 


18.48 


0.14 


5/7/5/7 


0.076 


0724+396 


2.753 


35.18 


19.15 


0.12 


19.93 


0.15 


18.78 


0.04 


18.75 


0.05 


12/12/12/12 


0.062 


0726+431 


1.072 


35.00 


19.52 


0.29 


20.25 


0.20 


18.72 


0.13 


- 


- 


6/8/6/- 


0.056 


0739+397B 


1.700 


35.11 


17.98 


0.14 


19.01 


0.07 


18.87 


0.14 


18.33 


0.05 


6/11/6/11 


0.050 


0740+380C 


1.063 


35.66 


17.59 


0.17 


18.23 


0.05 


17.65 


0.08 


17.51 


0.04 


5/7/5/7 


0.043 


0756+406 


2.016 


35.08 


19.08 


0.20 


19.80 


0.09 


19.24 


0.06 


18.64 


0.06 


12/12/12/12 


0.048 


0802+398 


1.800 


34.96 


19.26 


0.12 


20.01 


0.09 


19.97 


0.11 


19.70 


0.15 


12/12/12/12 


0.030 


0821+394 


1.216 


35.23 


17.51 


0.18 


18.21 


0.08 


17.39 


0.06 


17.25 


0.04 


5/7/5/7 


0.034 


0821+447 


0.893 


34.99 


17.42 


0.16 


18.28 


0.09 


17.64 


0.08 


17.56 


0.06 


5/7/5/7 


0.039 


0827+378 


0.914 


35.34 


17.78 


0.05 


18.37 


0.05 


17.94 


0.06 


17.71 


0.03 


5/11/5/11 


0.028 


0829+425 


1.056 


34.60 


18.69 


0.20 


19.17 


0.08 


18.93 


0.15 


18.47 


0.05 


6/11/6/11 


0.030 


0836+426 


0.595 


34.26 


20.63 


0.30* 


20.02 


0.17 


19.70 


0.30* 


19.23 


0.13 


5/7/5/7 


0.033 


0849+424 


0.978 


34.84 


18.59 


0.20 


19.72 


0.16 


18.73 


0.19 


18.85 


0.12 


6/8/6/8 


0.018 


0859+470 


1.462 


35.28 


17.93 


0.30* 


18.12 


0.04 


19.49 


0.30* 


18.07 


0.06 


5/7/5/7 


0.009 


0904+386 


1.730 


34.86 


17.73 


0.11 


18.37 


0.06 


18.11 


0.04 


17.77 


0.03 


12/12/12/12 


0.002 


0906+430 


0.670 


35.44 


19.18 


0.30* 


18.96 


0.09 


20.00 


0.30* 


18.22 


0.07 


5/7/5/7 


0.009 


0907+381 


2.16 


34.84 


16.41 


0.07 


17.58 


0.04 


17.47 


0.03 


17.22 


0.03 


11/11/12/11 


0.010 


0910+392 


0.638 


33.97 


18.42 


0.08 


19.22 


0.05 


18.81 


0.06 


18.58 


0.05 


12/12/12/12 


0.002 


0913+391 


1.250 


34.94 


18.70 


0.19 


19.54 


0.15 


18.80 


0.13 


18.84 


0.12 


6/7/6/7 


0.000 


0918+381 


1.108 


35.28 


20.06 


0.49 


19.81 


0.14 


19.23 


0.19 


18.62 


0.08 


6/7/6/7 


0.002 


0922+407 


1.876 


34.75 


19.52 


0.14 


20.59 


0.09 


20.42 


0.10 


19.76 


0.11 


11/11/11/11 


0.002 


0922+422 


1.750 


35.39 


18.33 


0.12 


18.92 


0.08 


18.36 


0.06 


17.96 


0.05 


7/7/7/7 


0.006 


0922+425 


1.879 


35.54 


19.94 


0.29 


21.12 


0.19 


21.01 


0.43 


18.95 


0.07 


6/7/6/7 


0.006 


0923+392 


0.698 


34.74 


15.72 


0.09 


16.89 


0.03 


16.62 


0.04 


16.28 


0.02 


6/7/6/7 


0.005 


0926+388 


1.630 


34.98 


19.12 


0.18 


19.53 


0.08 


19.42 


0.09 


19.29 


0.09 


12/12/12/12 


0.000 


0935+397 


2.493 


34.96 


21.15 


0.34 


21.34 


0.13 


20.89 


0.14 


20.52 


0.21 


12/12/12/12 


0.005 


0937+391 


0.618 


34.56 


17.60 


0.16 


18.28 


0.03 


18.23 


0.13 


18.04 


0.05 


6/8/6/8 


0.003 


0945+408 


1.252 


35.20 


17.00 


0.14 


18.05 


0.03 


17.70 


0.08 


17.58 


0.04 


6/8/6/8 


0.000 


0951+408 


0.783 


34.51 


18.76 


0.11 


19.56 


0.10 


18.77 


0.08 


18.92 


0.05 


7/11/7/11 


0.001 


0953+398 


1.179 


34.43 


18.59 


0.08 


19.56 


0.06 


19.14 


0.06 


18.73 


0.04 


10/10/10/10 
12/12/12/ 12 


0.000 


0955+387 


1.405 


35.29 


20.04 


0.33 


21.43 


0.41 


21.66 


0.89 


21.36 


0.34 


6/8/6/8 


0.004 


1007+417 


0.613 


34.87 


15.54 


0.07 


16.35 


0.02 


16.05 


0.03 


16.15 


0.02 


6/8/6/8 


0.000 


1015+383 


0.380 


33.65 


16.99 


0.06 


17.70 


0.03 


17.59 


0.05 


17.39 


0.05 


7/8/7/8 


0.000 


1020+400 


1.250 


35.01 


16.83 


0.11 


18.24 


0.08 


17.83 


0.09 


17.66 


0.06 


6/8/6/8 


0.000 


1030+415 


1.120 


34.70 


18.54 


0.11 


19.20 


0.08 


18.85 


0.14 


18.32 


0.05 


7/11/7/11 


0.000 


1105+392 


0.781 


34.87 


19.19 


0.34 


19.91 


0.08 


19.17 


0.20 


18.87 


0.07 


6/8/6/8 


0.000 


1109+437 


1.680 


35.99 


18.15 


0.14 


18.94 


0.06 


18.85 


0.13 


18.42 


0.05 


6/11/6/11 


0.000 


1111+408 


0.734 


35.55 


16.51 


0.08 


17.39 


0.03 


17.26 


0.06 






6/8/6/- 


0.020 


1116+392 


0.733 


34.01 


18.01 


0.12 


18.62 


0.06 


18.27 


0.04 


18.09 


0.04 


12/12/12/12 


0.000 


1123+395 


1.470 


34.63 


17.42 


0.09 


18.47 


0.05 


18.33 


0.04 


17.95 


0.03 


12/12/12/12 


0.000 


1128+385 


1.735 


34.44 


18.78 


0.23 


19.13 


0.06 


18.77 


0.06 


18.31 


0.03 


12/12/12/12 


0.000 


1141+400 


0.907 


34.28 


20.04 


0.20 


20.25 


0.07 


19.71 


0.08 


19.52 


0.08 


12/12/12/12 


0.000 


1142+392 


2.276 


35.33 


18.32 


0.15 


19.04 


0.06 


18.81 


0.06 


18.68 


0.05 


12/12/12/12 


0.008 


1144+402 


1.010 


34.43 


17.85 


0.19 


18.82 


0.04 


18.24 


0.08 


17.94 


0.04 


6/8/6/8 


0.000 


1148+387 


1.303 


35.27 


16.06 


0.10 


17.33 


0.04 


17.08 


0.05 


16.75 


0.03 


6/8/6/8 


0.000 


1148+477 


0.867 


35.05 


16.82 


0.08 


17.53 


0.03 


17.13 


0.03 


16.95 


0.02 


7/11/7/11 


0.004 


1203+384 


0.838 


34.42 


18.38 


0.10 






18.24 


0.05 


18.14 


0.04 


7/-/7/11 


0.000 


1204+399 


1.530 


34.47 


17.26 


0.08 


18.23 


0.05 


18.14 


0.04 


17.83 


0.04 


12/12/12/12 


0.009 


1206+439B 


1.400 


35.79 


17.77 


0.09 


18.45 


0.04 


17.97 


0.04 


17.34 


0.02 


7/11/7/11 


0.000 
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Table 2 — continued 



B3 name 


z 


log P 4 08 


U 


a(U) 


B 


a(B) 


V 


a(V) 


R 


a(R) 


Obs. date 


E(B - V) 


1228+397 


2.217 


35.26 


17.88 


0.08 


18.54 


0.06 


18.42 


0.06 


18.18 


0.07 


7/7/7/7 


0.000 


1229+405 


0.649 


34.31 


19.04 


0.12 


19.39 


0.10 


19.33 


0.14 


19.04 


0.07 


7/11/7/11 

1 / / / 


0.000 


1239+442B 


0.610 


34.37 


17.94 


0.07 


18.63 


0.03 


18.39 


0.04 


18.14 


0.03 


7,11/11/7,11/11 


0.000 


1240+381 


1.316 


34.29 


17.86 


0.06 


19.28 


0.05 


18.87 


0.06 


18.43 


0.05 


12/12/12/12 


0.000 


1242+410 


0.811 


34.68 


19.71 


0.20 






20.02 


0.15 


19.65 


0.13 


6/-/6/8 


0.000 


1247+450A 


0.799 


34.71 


17.05 


0.07 


17.93 


0.03 


17.62 


0.05 


17.47 


0.02 


7/11/7/11 


0.004 


1256+392 


0.978 


34.66 


19.22 


0.05 


19.62 


0.10 


19.12 


0.13 


18.85 


0.05 


7/11/7/11 


0.000 


1258+404 


1.656 


35.93 


18.32 


0.09 


18.63 


0.05 


18.45 


0.06 


18.18 


0.04 


7/11/7/11 


0.000 


1312+393 


1.570 


35.04 


21.41 


0.75 


22.22 


0.30 


21.97 


0.59 


21.76 


0.49 


11/11/6/11 


0.000 


1315+396 


1.560 


35.06 


17.75 


0.15 


18.69 


0.04 


18.61 


0.05 


18.20 


0.04 


11/11/11/11 


0.000 


1317+380 


0.835 


34.44 


19.65 


0.22 


20.42 


0.17 


18.98 


0.10 


19.64 


0.09 


7/11/7/11 


0.000 


1339+472 


0.502 


34.47 


20.13 


0.49 


20.26 


0.11 


19.23 


0.06 


18.74 


0.05 


6/11/6,12/11 


0.006 


1341+392 


0.768 


34.61 


20.68 


0.29 


21.18 


0.09 


20.80 


0.14 


20.23 


0.14 


11/11/7,11,12/11 


0.000 


1342+389A 


1.533 


35.15 


17.84 


0.11 


18.41 


0.04 


18.32 


0.04 


17.97 


0.03 


11/11/11/11 


0.000 


1343+386 


1.844 


35.31 


17.55 


0.10 


18.10 


0.04 


17.87 


0.03 


17.50 


0.03 


11/11/11/11 


0.000 


1348+392 


1.580 


34.88 


19.19 


0.14 


19.99 


0.09 


19.79 


0.08 


19.61 


0.09 


12/12/12/12 


0.000 


1355+380 


1.561 


34.93 


18.65 


0.07 


19.48 


0.04 


19.41 


0.07 


18.78 


0.05 


11/11/11/11 


0.000 


1357+394B 


0.804 


34.20 


18.08 


0.06 


18.85 


0.05 


18.53 


0.05 


18.31 


0.04 


12/12/12/12 


0.000 


1416+400 


0.473 


33.88 


19.89 


0.15 


20.37 


0.05 


19.94 


0.11 


19.40 


0.09 


11/11/11/11 


0.000 


1417+385 


1.832 


34.41 


18.54 


0.07 


19.19 


0.06 


19.03 


0.05 


18.65 


0.06 


12/12/12/12 


0.000 


1419+399 


0.622 


33.35 


18.67 


0.12 


19.00 


0.06 


19.01 


0.07 


18.75 


0.03 


12/12/12/12 


0.000 


1435+383 


1.600 


34.89 


17.47 


0.09 


18.20 


0.05 


18.04 


0.04 


17.78 


0.04 


12/12/12/12 


0.000 


1444+417A 


0.675 


34.60 


18.73 


0.39 


19.42 


0.21 


19.39 


0.34 


19.08 


0.12 


6/11/6/11 


0.002 



* indicates unreliable photometry (second part of the night of Feb 5th) 
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Figure 1. Individual rest- frame SEDs of the B3-VLA quasars plotted in order of increasing redshift. Vertical dashed lines indicate typical 
broad emission lines observed in quasars: H/3, MgoA2798, Coi]A1909, ClvA1549 and Lya. 
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Figure 1 - continued 
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Figure 2. Histograms of the UBVR magnitudes for the 64 
quasars with good photometry 



from the literature, which can be used to estimate their av- 
erage contribution to the broad-band fluxes. The listed EWs 
include the values obtained by Baker & Hunstead (1995) for 
Molonglo radio quasars, which is a sample selected at 408 
MHz, with 5*408 > 0.95 Jy, and therefore appropiate for com- 
parison with our sample (although it does not list the EWs 
for the Fell bumps). Also listed are the EW measurements 
for radio- loud quasars given by Zheng et al. (1997, from an 
UV composite which extends to 3000 A) and the EWs for 
optically-selected LBQS quasars from Francis et al. (1991) 



and P.J. Green (1998). Francis et al. (1991) mention that the 
height of the bumps derived from their work would be de- 
creased if the regions around 1700 A and 2650 A were taken 
as continuum, although at the cost of a strongly curved con- 
tinuum. Green (1998) uses as continuum the regions 2645- 
2700 A, 3020-3100 A and other wavelengths around 4430 
and 4750 A, and finds that the Feii+Balmer bump in the 
range 3100-3800 is absent in his composite. The EW for Fe 
II 2400 from Green (1998) is very similar to that measured 
by Zheng et al. (1997). 

Baker & Hunstead (1995) obtain qualitative estimates 
of the contribution of Fell bumps for different radio mor- 
phologies from the height of the composite spectra relative 
to a power-law continuum from ~2000 to ~5000 A. The 
authors conclude that the Fell bumps are absent in lobe- 
dominated and compact-steep-spectrum quasars, but rather 
strong in core-dominated quasars, which resemble in vari- 
ous line properties optically-selected quasars. The contin- 
uumm of optically selected quasars is known to curve at 
around 3000 A (Francis et al. 1991, Natali et al. 1998), and 
the strength of the bumps obtained by Baker & Hunstead 
will be significantly decreased if the underlying continuum 
was allowed to curve. Boroson & R.F. Green (1992) had 
previously found that flat-spectrum quasars have stronger 
Fen 4434-4684 A than steep-spectrum ones, and that, as a 
whole, radio-loud quasars have fainter emission than radio- 
quiet ones. The latter result was also reported by Cristiani & 
Vio (1990); Fell 2400 is the only Fell bump revealed in their 
composite spectra for radio-loud quasars, and the feature, 
almost absent, is less pronounced than in the composite for 
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Table 3. Average equivalent widths of quasar emission lines 





Baker &; 
Hunstcad 


2 lie ii u 
ct al. 




Fr£Liicis 
ct al. 


Adopted 


Civ 


93 


77 


68 


37 


93 


Cm] 


23 


17 




20 


23 


Mgn 


58 


50 


53 


50 


58 


Fen 2400 




38 


39 


35 


38 


Fen 3400* 






absent 


76 


38 


Fen 4500 






29 


41 





* Includes Balmer emission 

radio-quiet quasars. Zheng et al. (1997) derive however a 
lower EW for Fen 2400 for radio-quiet quasars (22 A versus 
38 A), showing evidence of the uncertainties in the contri- 
bution of these bumps. 

Adopting as EWs for the Civ, Cm] and Mgn lines the 
values from the Molonglo sample we infer the maximum con- 
tribution to the broad-band fluxes for the Civ line, amount- 
ing on average to 25 per cent if the line is included in the 
U, B or V band. For the Cm] line the average contribu- 
tion is below 7 per cent and for Mgn in the range 6-11 per 
cent. As EW for Fell ~ 2400 A we adopted the value for 
the radio quasar sample by Zheng et al. (1997). The EWs 
for Civ, Cm] and Mgn by Zheng et al. (1997) are in fact 
very similar to those from the Molonglo quasars, and the 
radio sample is also very similar to ours in terms of opti- 
cal luminosity. From the average EW of 38 A the inferred 
contribution to the broad-band fluxes is 5-9 per cent. The 
Feii+Balmer ~ 3400 feature lies outside the spectral region 
covered by the composite by Zheng et al., and is absent in 
the composites by Green (1998) and Cristiani & Vio (1990). 
Adopting half the EW from Francis et al. (1991), i.e. around 
38 A, its expected contribution would be below 6 per cent 
when included in the V or R bands. 




Figure 3. Distribution of the broad-band colours U — R, B — R 
and U — V as a function of redshift. The underlined symbols 
correspond to the red quasars B3 0918+381 and B3 1339+472 
(see Section 6.1) 



4.2 Optical broad-band colours of the quasars 

The distribution of the U — R, B — R and U — V colours 
as a function of redshift is shown on Figure 3. The broad- 
band colours do not show a clear variation with redshift. 
If any, this would be a blueing of the U — V colour with 
redshift up to z around 2.5. The mean colours derived for 
these quasars (indicated on the figure with dotted lines) are 
(U-R) = -0.08, (B~R) — 0.64 and (U — V) = -0.38 with 
dispersions 0.45, 0.27, and 0.42, respectively, and correspond 
to observed spectral indices a obB tj-r = —0.50, a obs b-r = 
—0.39 and a obs u-v = —0.75 with dispersions 0.64, 0.53, 
and 0.96. From the comparison of the B — R and U—V 
spectral indices a trend is found in the sense that the SED is 
steeper at higher frequencies, in agreement with the results 
found by Natali et al. (1998). 



5 SPECTRAL ENERGY DISTRIBUTION OF 
THE INDIVIDUAL SOURCES 

In this section we discuss the shape of the individual SEDs 
of the B3 quasars. The SEDs of the 61 sources with available 
photometry at the four bands were fitted through \ 2 mini- 
mization, using a power-law model and a quadratic model. 
The second model was chosen as a simple representation for 



the SEDs curved on the logSV-logf plane. A fit was accepted 
if the probability Q that the \ 2 should exceed the particu- 
lar measured value by chance was higher than 1 per cent. 
The best-fit spectral indices and their errors are listed in Ta- 
ble 4. Forty-one quasars have acceptable fits as power-laws. 
Thirty-nine of these have also acceptable fits as quadrat- 
ics and in 6 cases the quadratic model gives a significant 
improvement of the fits (higher than 85 per cent using an 
F-test). These six cases are labelled as p^q on the last col- 
umn of Table 4. There are cases where the shape of the 
SED is clearly curved, but a power-law model is acceptable 
due to large errors. Similarly, SEDs which resemble power- 
laws to the eye and have small photometric errors may not 
have acceptable fits for this model. The two sources with 
acceptable fits as power-laws but not as quadratics are B3 
0849+424 and B3 1144+402, and both sources have poor 
power- law fits (Q < 0.02). The SEDs of these sources could 
be contaminated by emission lines at the U and V bands. 

Ten additional quasars without acceptable fits as power 
laws can be fitted with quadratics. The remaining 10 sources 
do not have acceptable fits with any of the models (indi- 
cated with an hyphen in the last column of Table 4). For 
some of these sources the lack of good fits could be due 
to the presence of emission peaks in their SEDs related to 
contamination by emission lines. We note for instance that 
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Table 4. Results of power-law and quadratic fits 



B3 name 


a 


a(a) 


Model 


B3 name 




a(a) 


Model 


B3 name 


a 


a(a) 


Model 


0701 -4-^Q9 


0.22 


0.25 


1 


nQ354-3Q7 


-1.10 


0.41 


P 


1 9984-3Q7 


-0.10 


0.13 


P 


07044-^84 


-0.81 


0.25 


P 


0Q374-3Q1 


0.37 


0.11 


P 


1 9904-405 


-0.32 


0.17 


q 


07944-^Qfi 


AQ 


1 fi 
U. 1U 






04 


OR 


P 


1 93Q4-449R 


17 

-U. J. i 


OS 

u.uo 


p— >q 


073Q4-^Q7R 


-0.18 


0.14 


p— >q 


U3J I TWO 


-0.27 


0.14 






-0.09 


0.10 




07404-^SOr 1 


-0 41 


12 


q 




^fi 

-u.ou 


10 
U. 1U 


P 


1 9474-450 A 


008 

-u.uuo 


07 

u.u / 


p 


07^fi4-40fi 


1 1C 
-l.lo 


U.lO 


p 


1 0074-41 7 


44 


O^i 

u.uo 




1 95fi4-3Q9 


111 


10 
u. 1U 


q 


08094-^QR 


0.19 


0.25 


p 


iuiotooo 


0.11 


0.09 


P 


1 9584-404 


-0.29 


0.10 


q 


0821 4-^Q4 


-0.68 


0.15 




1 020+400 


0.20 


0.15 






-0.09 


0.11 


p 


0R91 -4-447 


-0.28 


0.16 




1 0^04-41 ^ 


-0.90 


0.13 


P 


1 ^^Q4_479 


-1.73 


0.20 


q 


0R974-^7R 


-0.52 


0.07 


q 


1 1 0^4-3Q9 


-1.16 


0.20 


P 


1 ^41 4-^Q9 


-1.06 


0.29 


p 


08904-49^ 


-0.53 


0.17 


p 


1 1 004-4^7 


-0.18 


0.14 


P 


1 ^494-^SQ A 


-0.13 


0.09 


p 


0S4Q4-494 


-0 41 


9Q 


p 




-0 22 


11 

U. 1 J. 


q 


1 ^4^4-^Rfi 
IOIOtOOU 


4^ 


OS 

u.uo 


p 


0904+386 


-0.37 


O.lf 


p 


1123+395 


-0.07 


0.09 


p^q 


1348+392 


0.02 


0.19 


p 


0907+381 


0.36 


0.08 


p 


1128+385 


-0.80 


0.12 


p 


1355+380 


-0.41 


0.11 


p 


0910+392 


-0.40 


O.lf 


p 


1141+400 


-0.83 


0.19 


q 


1357+394B 


-0.27 


0.10 


p 


0913+391 


-0.52 


0.28 


p 


1142+392 


0.12 


0.14 


p^q 


1416+400 


-1.14 


0.19 


p 


0918+381 


-1.61 


0.28 


p 


1144+402 


-0.85 


0.11 


p 


1417+385 


-0.36 


0.12 


p 


0922+407 


-0.38 


0.22 


p^q 


1148+387 


-0.14 


0.08 




1419+399 


0.08 


0.11 




0922+422 


-1.02 


0.14 


P 


1148+477 


-0.25 


0.06 


q 


1435+383 


-0.05 


0.01 


p 


0923+392 


-0.22 


0.06 




1204+399 


0.15 


0.10 


p^q 










0926+388 


0.12 


0.20 


P 


1206+439B 


-1.28 


0.08 


p 











It 



15 



10 



Figure 4. Spectral index versus rcdshift for the 6f quasars with 
good photometry and available data in the four bands. Filled 
circles correspond to power-law fits formally acceptable 



-1.5 



-1 -0.5 
spectral index 



0.5 



Figure 5. Spectral index distribution for the quasars shown on 
Fig. 4. The dotted curve shows the gaussian fit for the quasars 
with a > -0.9 



contamination by Civ could be related to the maxima in the 
U band for B3 1020+400, B3 1148+387 and B3 1240+381, 
and in the V band for B3 0724+396. 

Figure 4 shows the distribution of spectral indices ver- 
sus redshift for the 61 sources. The slopes were obtained 
from fixed observed wavelengths, therefore they correspond 
to different rest-frames (from 2600-5000 for z = 0.4 to 1000- 
1900 for z — 2.8). The slopes show a wide range, with values 
from 0.4 to -1.7, regardless that the fits are formally accept- 
able or not. The mean and the dispersion for the total sam- 
ple are -0.39 and 0.38 respectively (-0.41 and 0.40 for the 
sources with formally acceptable power-law fits). The com- 
parison of these values with those obtained from two-band 
colours (Sect. 4.2) show the best agreement for B — R, with 



Qobs b-r = —0.39 and standard deviation 0.53. The reason 
for this is that the B and i?-band errors are low, compared 
to those for the [/-band, and the errors are weighted for the 
power-law fits. 

Although a mean spectral index was obtained from the 
sample, the distribution of slopes is asymmetric, showing a 
tail to steep indices (see Fig. 4 and the histogram on Figure 
5). Considering the total sample of 61 sources, the distribu- 
tion of spectral indices for a > —0.9 is well represented by 
a gaussian with a mean of -0.21 and a dispersion of 0.34. 



6 SPECTRAL ENERGY DISTRIBUTION 
FROM THE COMPOSITE SEDS 
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6.1 Normalized composite SEDs 

The shape of the SED of the quasars can be analysed 
through a "composite spectrum", in which the individual 
SEDs are merged in the rest-frame, adopting a specific cri- 
terium for the scaling of the fluxes. The coverage of the 
optical photometry and the range of redshifts of the quasars 
allows to study their spectrum in the range ~1300 — 4500 
A. We chose to normalize the individual SEDs to have the 
same flux density at a fixed wavelength A n within the ob- 
served range, and the flux density at A n was obtained by 
linear interpolation between the two nearest data points. 
With this procedure we could not have an appropriate nor- 
malization for the whole sample (the coverage for the lowest 
redshift quasars is 2500 - 5000 A and for the highest red- 
shift ones 1000 — 2000 A), therefore we obtained composite 
SEDs for different normalization wavelengths. The selected 
An were 3800, 3500, 3200, 2400, 2200 and 2000 A and the 
corresponding normalized SEDs are shown on Figure 6. The 
use of a broad range of values for A n is appropriate to anal- 
yse the possible dependence of the SED shape with A n . The 
low separation between the A n values, typically around 300 
A, allows for a large overlap between the composites and 
for their comparison in a continuous way. A normalization 
around 2800 A was not considered, to avoid the Mgn emis- 
sion line. The fluxes at the normalization A n =2000 could 
be contaminated by the weaker Cm]A1909 line, and those 
at the normalizations 2400, 3200 and 3500 by the emission 
bumps at 2250-2650 and 3100-3800. However, we have seen 
in Sect. 4 that the expected contributions of these features 
are weak; typically below 10 per cent. A 10 per cent con- 
tribution corresponds to a vertical shift in logS,, (Figure 6) 
of 0.04, which is clearly lower than the dispersion of the 
spectral points in the composites. 

Some of the spectral points clearly deviating in the 
composite SEDs correspond to the two reddest quasars, B3 
0918+381 and B3 1339+472 (a < -1.6). The extreme red 
colours of these two sources are uncommon in the sample, 
as can be seen from Fig. 3, where the symbols for the two 
sources appear underlined. Both sources show curved spec- 
tra in the logSV-logA plane. B3 0918+381 has an acceptable 
fit for the power-law model, but this is due to the large mag- 
nitude error at the [/-band. On Fig. 6 the normalized SEDs 
of these sources are indicated over the remaining spectral 
points, showing the notorious discrepancy relative to the 
average composite SEDs. These two sources with peculiar 
SEDs will not be considered for the discussion in this and 
the next section. 

The overall shape of the composite spectra is found to 
be very uniform for the different normalizations. The dis- 
persion of the spectral points is artificially reduced near 
the wavelength A n . The spectral points appear to trace pre- 
dominantly the continuum; only the CivA1549 line appears 
prominent on the composites, as an "emission feature". No 
other emission bumps/lines are apparent from the compos- 
ites, in agreement with the low contributions expected (Sect 
4.1). The Civ feature is revealed in the three composites 
covering its wavelength and arises from the [/-band data of 
various quasars, most of which were mentioned in Sect. 5 to 
have likely contamination by this line. The emission feature 
appears in the broad-band composite displaced to the red, 
peaking at around 1610 A. The displacement is around five 



times lower than the rest-frame width over which the line is 
detected. 

A general trend on Fig. 6, apparent especially from pan- 
els (b) to (d), is a steepening of the spectrum from large to 
short wavelengths, ocurring at around 3000 A. Above this 
wavelength the SED appears to be rather flat. This trend 
is less evident on panel (a), probably because of the small 
number of spectral points at A < 3000 A. The trend is weak 
in panel (e) and practically disappears in panel (f), but here 
the number of points with A > 3000 A is small. Panels (d) 
and (e) appear to show that the SED flattens again below 
2000 A. This trend practically dissapears in panel (f), in 
spite of the large number of points, thus this flattening is 
not so clear in principle as the steepening at 3000 A. 

6.2 Power-law fits 

We have obtained power-law fits of the composite SEDs in 
the regions above and below the 3000 A break. The selected 
ranges were 4500-3000 A (referred to as "blue") and 2600- 
1700 A (referred to as "UV"), which roughly correspond to 
the regions limited by H/3 and MgnA2798 and MgnA2798 
and CivA1549. These spectral points were taken as repre- 
senting the continuum, since the contamination by emission 
lines/bumps in this range is expected to be very weak. The 
results of the fits, obtained by least-squares minimization, 
are listed in Table 5, including for each composite the slopes 
and their errors, the number of spectral points used and the 
redshift range of the quasars. Slopes with an asterisk cor- 
respond to lower quality fits, since the spectral points do 
not cover the whole wavelength range. Figure 7 shows the 
spectral index values as a function of A n , illustrating its vari- 
ation between the different composites. The four composite 
SEDs for which power-laws were fitted in both ranges show 
a steepening from low to high frequency. The differences 
Qbiuo — «uv for these composites and their errors are listed 
in Table 5. For the following discussion only the good quality 
fits will be considered. 

Concerning the 4500—3000 A range, the spectral in- 
dices for all the normalizations are consistent within their 
errors. The first three normalizations, with A n =3800, 3500 
and 3200, show a better agreement with each other. Aver- 
aging the spectral indices for the fits for A n =3800, 3200 and 
2400 (the fit for 3500 A has a large overlap with those for 
3200 and 3800 A) we obtain (a b i U o) = 0.11 ± 0.16. 

The high frequency fits have again best-fit slopes con- 
sistent with each other within the errors. For the high fre- 
quency normalizations at A n =2400, 2200 and 2000, the spec- 
tral indices are very similar, although the range of spec- 
tral points used do not show a wide overlap. The spectral 
index for the low frequency normalization at A n =3200 is 
steeper. Averaging the slopes obtained for A n =3200, 2400, 
2200 and 2000 we find (aw) = -0.66 ± 0.15. Consid- 
ering the average measured values of Obiue and aw and 
their errors we find a steepening towards high frequencies 
(ablue) - (auv) = 0.77 + 0.22. 

Although the expected contribution of emission 
lines/bumps in the ranges used for the power law fits is 
low, it is interesting to analyse the possibility that the slope 
change is artificially produced by contamination due to the 
Fen bumps at 3100-3400 and 2250-2650, since these fea- 
tures would enhance the emission used for the power-law 
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Figure 6. Composite SED for different normalization wavelengths. Spectral points from different optical bands are shown with different 
symbols. Filled squares, empty squares, filled triangles and empty triangles correspond respectively to R,V,B and [/-band data. The 
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Table 5. Power-law fits in the intervals 4500-3000 A and 2600-1700 A for quasars in the redshift range 0.38-2.49 



An 


"blue 


IT 


N 


z range 


"uv 


a 


N 


z range 


O^blue — "UV 


a 


3800 


0.20 


0.12 


34 


0.38-0.84 














3500 


0.23 


0.11 


41 


0.38-1.01 


-0.73* 


0.35 


38 


0.47-1.01 


0.96* 


0.37 


3200 


0.21 


0.12 


45 


0.38-1.18 


-0.82 


0.25 


47 


0.47-1.18 


1.03 


0.28 


2400 


-0.08 


0.23 


50 


0.57-1.32 


-0.61 


0.11 


96 


0.57-1.88 


0.53 


0.25 


2200 


-0.25* 


0.35 


36 


0.69-1.32 


-0.57 


0.10 


93 


0.69-2.16 


0.32* 


0.36 


2000 










-0.66 


0.10 


80 


0.83-2.49 







fits (4500-3000 and 2600-1700 A) in the regions next to 
the break, producing a slope change in the observed sense. 
However, a power law fit for the whole range from 1700 
to 4500 A for the composites (b) to (d), where the break 
is more obviously detected, yields a clear excess emission 
only in the region 2650-3200 A i.e. between the two bumps. 
In fact some excess there could be due to Mgn, but if the 
Fen bumps were responsible for the slope change, the ex- 
cess emission should extend over the Fen ranges 2250-2650 



and 3100-3400. Therefore, although some of the broad-band 
spectral points are expected to include contamination due to 
emission lines and/or bumps, the break in the overall SED 
detected at ~ 3000 A is most likely related to an intrinsic 
change in the continuum of the quasars. A steepening of 
the quasar's continuum at 3000 A was reported by Natali et 
al. (1998) from spectra of optically selected quasars, and is 
also found in the composite spectrum of optically selected 
quasars by Francis et al. (1991). The slopes obtained by Na- 
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tali et al. (1998) were (a Wuc ) ^ 0.15 for the range 2950-5500 
A, and (auv) — ~0.65 for the range 1400-3200 A, in good 
agreement with our values. 



6.3 Redshift dependence of the blue/UV 
continuum shape 

A striking characteristic of the best-fit slopes on Table 5 and 
Fig. 7 (top panels) is how they keep roughly constant for 
the normalizations at either side of the break around 3000 
A, suffering the largest variation when the normalization 
moves from one side of the break to the other. In particu- 
lar, the change in auv from A n =3200 to lower values of A n 
is consistent with the description outlined on Section 6.1 of 
a possible flattening at around 2000 A for the composites 
with A n =2400 and 2200 A. It is important to note that the 
change in normalizations from A n > 3200 A to A n < 2400 A 
in our study implies a substantial change in the redshifts of 
the quasars whose spectral points enter the fits. For instance, 
whereas that the high frequency fit for A n =3200 comprises 
47 spectral points with 0.47 < z < 1.18, the next compos- 
ite, with A n =2400, includes all but one of these points plus 
50 additional spectral points from 1.21 < z < 1.88. For the 
final composite, with A n =2000, only 23 out of the 80 spec- 
tral points used for the fit are in common with the fit for 
A n =3200. A similar variation of the spectral points/redshifts 
with An occurs for the low frequency fits, although in this 
case the low redshift limit for A n =3800, 3500 and 3200 is 
similar and the variations of spectral points/redshifts arc 
lower for all the normalizations (see Table 5) . It is therefore 
interesting to analyse whether a dependence of the shape of 
the continuum spectrum with redshift is present. 

From the results of the fits (Table 5, Fig. 7, and see 
also Fig. 6) we decided to separate the quasar sample in 



a low redshift bin with z < 1.2 and a high redshift bin 
with z > 1.2, and perform the power-law fitting for the two 
subsamples. With the selected limit, all the spectral points 
used for the normalizations with A n > 3000 A correspond 
to the low redshift bin. The composites with A n < 2400 A 
include spectral points from quasars in the two redshift bins. 
Figures 8 and 9 show the composite SEDs for the two quasar 
subsamples. The results of the power-law fits are given in 
Table 6 and the variation of the slopes as a function of A n 
is shown on Fig. 7. 

Concerning the low redshift quasars, we find now a 
better agreement between the slopes of the low-frequency 
fits obtained for the normalizations on either side of the 
break. Averaging the values for the normalizations with 
A n = 3800, 3200 and 2400 we find (a b i U e) = 0.21 ± 0.16. 
For the high frequency fits there is also a better agree- 
ment in Quv between the normalizations at either side of 
the break. The agreement is due to the steepening of auv 
for the normalizations with A n < 2400 A when high redshift 
quasars have been excluded. The average value of auv is 
now (auv) = — 0.87T0.20. We find now a larger change from 
Qbiuc to auv, with (abi uo ) — (auv) = 1.08 ±0.26. The larger 
steepening is evident from the comparison of the top and 
bottom panels of Fig. 7. The good agreement in the slopes 
obtained for the different normalizations indicates that the 
selection of A n does not affect the derived slopes. 

For the high redshift quasars only the high frequency re- 
gion can be studied (see Fig. 9). The average spectral index 
from the three normalizations is (auv) = — 0.48T0.12. Com- 
paring the values of auv for low and high redshift quasars 
we find a difference A(auv) = 0.39 ± 0.23. This difference 
although weak, is larger than the error, and is clearly evident 
from the results on Table 6 and Fig. 7, and the comparison 
of figs. 8 and 9. The flattening below around 2000 A noted 
on Sect. 6.1 for the high frequency normalizations is related 
to this hardening of auv for high redshift quasars. 

Our result of a dependence of the spectral index with 
redshift, in the sense that quasars at z > 1.2 tend to be flat- 
ter, i.e. harder, than low redshift ones, in the range 2600- 
1700 A, is similar to that found by O'Brien et al. (1988) 
for an IU E-selected sample, with slopes in the region 1900- 
1215 A ranging from -0.87 for z < 1.1 to -0.46 for z > 1.1. 
We note the remarkable agreement between these values and 
those obtained for our sample. This trend with redshift was 
not confirmed by Natali et al. (1998) for their sample of 
optically selected quasars, although their limits for the de- 
tection of a variation were around 0.4, which is comparable 
to the measured change in the slope found in O'Brien et al. 
and in this work. 

A possible concern on our results is whether the 20 
per cent incompleteness of the sample could produce a bias 
against the inclusion of red quasars at high redshifts. How- 
ever, the distribution of R magnitudes in the range 2=0.5- 
2.5 is roughly constant, and does not suggest that the frac- 
tion of missed quasars (R > 20) should be higher at the 
higher redshifts. On the other hand, 5 good B3-VLA quasar 
candidates excluded from the B3-VLA Quasar Sample were 
present on the blue POSS-I plate but not on the red POSS 
plate. Since the magnitude limit for POSS blue is around 
20.5-21, these quasars have B — R colours similar or bluer 
than the average value for the quasars in the sample (B — R— 
0.64, Sect. 4.2). Although this number of quasars is small, 
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Table 6. Power-law fits for quasars in the redshift range 0.38-1.18 (upper part) and 1.21-2.49 (bottom part) 
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3500 
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it illustrates the presence of blue or normal colours among 
the missed quasars. Since the incompleteness does not ap- 
pear to be biased towards the exclusion of red quasars at 
high redshifts, the «uv — z trend is likely a real effect, not 
originated by incompleteness. 



6.4 Relation of the blue/UV continuum shape to 
radio power 

Figure 10 presents the P408 — z diagram for the 73 quasars 
observed for this work. In agreement with the expectations 
for flux-limited samples, high redshift quasars tend to have 
higher radio powers than low redshift ones, although the B3- 
VLA Quasar Sample is not strictly a flux-limited one (see 
Sect. 2). Using Spearman's correlation coefficient r s , we find 
for the trend r s =0.52, with a significance level P > 99.999. 
In order to check whether the Quv — z trend could arise 
from an intrinsic auv — P408 correlation and the Ptos — z 
trend, we performed an analysis similar to the one presented 
in sections 6.2 and 6.3, obtaining separated power-law fits 
for a "high-radio-power subsample" and a "low-radio-power 
subsample". Instead of considering several composites, the 
dependence on radio power was analysed using only the com- 
posite with A n =2400. This composite includes the largest 
number of quasars, and the redshift distribution comprises 
low as well as high redshift ones, with a median z of 1.15, 
similar to that of the whole sample. 

The composite with A n =2400 includes 54 quasars (ex- 
cluding B3 0918+381). The median of P408 is 10 34 ' 85 erg s" 1 
Hz -1 and this value was used as the division limit between 
low and high radio power. The low-power quasars have a 
mean redshift of 1.00 and a mean radio power IO 3448 erg s _1 
Hz -1 , and the same parameters for the high-power quasars 
are 1.35 and 10 35 19 erg s _1 Hz -1 . The results of the derived 
spectral indices Quv for the two subsamples are listed on 
Table 7. The table includes for comparison the Quv values 
obtained for all the quasars and for the subsamples sepa- 
rated by redshift, for the same composite. The third column 
of the table gives Aquv and its error for each couple of sub- 
samples with a varying parameter (z or Pios), and the last 
column gives the probability that the difference in slope is 
statistically significant. 

The largest difference in auv and the highest signif- 
icance is found for the subsamples separated by 2, with 
Aauv = 0.54 ± 0.22, and 98.84 per cent significance (F- 
test). For the subsamples separated by P408 the difference 
in auv is lower than the errors. The trend auv — z is there- 



36 



so 35 



• • • * 



34 - 



33 



Figure 10. -P408 MHz versus redshift for the 73 B3-VLA quasars 
observed for this work 



Table 7. ojtjv f° r the composite with A n =2400 and different sub- 
samples, separated by redshift, radio power and blue/UV lumi- 
nosity 





«uv 


A(q uv ) 


Significance of 
the slope change 


All 2 


-0.61 ±0.11 






z < 1.2 
z > 1.2 


-0.94 ±0.17 
-0.40 ±0.14 


0.54 ±0.22 


98.84% 


P408 < io 34 - 85 
P408 > io 34 - 85 


-0.71 ±0.18 
-0.51 ±0.14 


0.20 ±0.23 


27% 


L2400 < 10 30 - 67 

L2400 > io 30 - 67 


-0.76 ±0.18 
-0.46 ±0.13 


0.30 ±0.22 


89% 



fore unlikely to be a secondary correlation from an intrinsic 
Quv — P408 correlation in conjunction with the P 40 8 — 2 trend. 
Moreover, the last trend is less significant for the group of 54 
quasars in the A n =2400 composite (r a =0.44 and P=99.906). 
This weak trend explains in fact the reverse result, i.e. that 
the auv — 2 trend in combination with the P408 — 2 trend 
does not cause a auv — P408 trend. 
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Figure 11. L ° versus z for the 54 quasars with measured 
b 2400 A H 

(interpolated) flux density at 2400 A 

6.5 Relation of the blue/UV continuum shape to 
blue/UV luminosity 

The possible relation between auv and the blue/UV lumi- 
nosity was also investigated. In order to construct each com- 
posite SED we had to interpolate the flux density at the nor- 
malization wavelength A„, and therefore we could obtain the 
monocromatic luminosity at A n straightforward. The possi- 
ble trend Lbiue/uv — Q uv was analysed in the same way as 
for the P408 — auv trend, measuring auv in the composite 
with An = 2400 for two subsamples separated by L2400 • The 
separation value was set at L2400 = 10 30 ' 67 erg s" 1 Hz -1 , 
which is the median value for the quasars in the composite. 
The values of auv (Table 7) indicate a weak flattening with 
luminosity, only slighly larger than the errors, but signifi- 
cant from the F-test, with P=89 per cent. The significance 
is lower than that found for the auv — z trend. 

Figure 11 shows L2400 versus z for the quasars in the 
same composite. A clear correlation is found, in the sense 
that high redshift quasars are more luminous, similar to 
the trend with radio power. A Spearman correlation gives 
r s =0.52 and significance 99.994 per cent. Therefore the three 
parameters auv, z and L2400 are correlated with each other, 
being in principle impossible to determine which are primary 
and which are secondary correlations. The quasars with lu- 
minosities below the limit of 10 30 67 erg s" 1 Hz" 1 have a 
mean redshift and luminosity of 1.01 and 10 30 27 erg s" 1 
Hz" 1 , and the same parameters for the luminous quasars 
are 1.33 and 10 31 02 erg s" 1 Hz" 1 . 

Wandel & Petrosian (1988) obtained predictions of the 
expected UV slope versus UV luminosity for various accre- 
tion disc models, spanning a range of values for the black 
hole mass, accretion rate and viscosity parameter. They 
found that the effect of a flattening of the spectrum for the 
more luminous quasars (or at higher redshifts), reported by 
O'Brien et al. (1988), could be explained as evolution along 
curves of constant black hole mass, with accretion rate de- 
creasing with time. The emitted spectrum for these mod- 
els under the assumption of constant black hole mass is 
less luminous and more steep in the UV as the accretion 



rate decreases, and a decrease in accretion rate with time 
could explain the observed Luv — z and auv — z trends. 
Our data show the same trend of a flattening of auv with 
z and Luv, as well as a Luv — z correlation, therefore the 
same interpretation could apply. The diagrams auv — iuv 
presented by Wandel & Petrosian (1988) correspond to a 
slope in the range 1500 — 1000 A and the luminosity at 
1450 A. We have represented in these diagrams the average 
slope for low-redshift and high-redshift quasars (separation 
at z = 1.2), using the composite with A n = 2400, and the 
2600-1700 A slopes listed in Table 7. The luminosities were 
obtained transforming the average values at 2400 A to 1450 
A, adopting the measured slopes. Although we use a slope 
different than that of the models, at least for high redshift 
quasars the slope does not apperar to change between the 
two ranges (Fig. 9). A straight line connecting the location 
of the two subsamples on the auv — iuv plane would run 
approximately parallel to the curves of constant black hole 
mass, with m increasing in the same direction as the red- 
shift. Therefore our results appear to be consistent with the 
interpretation outlined above, and the comparison with the 
models yields a roughly constant black hole mass, in the 
range M ~ 1O 8 ' 2_8 ' 5 M0, and accretion rates (in units of the 
Eddington value) ranging from rh ~ 1 for the quasars with 
z = 0.6 — 1.2 to m ~ 10 for the quasars with z = 1.2 — 1.9. 



7 THE CORRELATION BETWEEN BLUE/UV 
LUMINOSITY AND RADIO POWER 

Figure 12(a) shows L2400 versus P408 for the 54 quasars of 
the A n =2400 composite. L2400 and P408 MHz are strongly 
correlated, with r s =0.62 and P > 99.9999 per cent. The op- 
tical completeness of the sample is high and we do not expect 
important selection effects in the optical that could explain 
the correlation. The correlation could be in principle artifi- 
cially induced by independent evolution of both parameters 
with redshift. In fact the quasars with high blue/UV lumi- 
nosity tend to have large redshifts and those weaker tend to 
have low redshifts, and a similar behaviour occurs for the 
radio power. However, this alternative is ruled out by the 
redshift distribution on Figure 12(a), showing that the half 
subsample with lower redshifts shows a similar L2400 — P408 
trend as the whole sample (with unrestricted 2). The cor- 
relation is confirmed in the flux-flux plane on Figure 12(b), 
with r a = 0.47 and P = 99.966 over two decades in both 
radio and blue/UV flux. 

The most likely origin for the radio-optical correlation 
is therefore a direct link between the blue/UV emission of 
the quasar nucleus, related to the accretion process, and the 
radio synchrotron emission at 408 MHz. A similar corre- 
lation was reported by Serjeant et al. (1998) on the basis 
of complete samples of radio-steep-spectrum quasars (from 
Molonglo, 3CR and the Bright Quasar Survey), in the red- 
shift range 0.3-3. A linear fit for our data in the logI/2400- 
logP 40 8 MHz plane yields L2400 oc P^mhz > witn a disper- 
sion in blue/UV magnitudes of ~ 0.9, for quasars in the red- 
shift range 0.6 — 1.9. The best-fitting slope is very similar to 
that obtained by Serjeant et al. (1998), with Lb oc ^los mhz> 
although these authors give a larger dispersion, ~1.6 mag. 

In the light of an intrinsic correlation between radio 
power and blue/UV luminosity we can investigate in more 
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Figure 12. Panel (a): L2400 versus P408 f° r 54 B3-VLA quasars in this work. Panel (b): S2400 versus S408 flux densities for the same 
quasars. For the radio fluxes we used the fc-correction C = — (1 + a) log(l + 2), where a is the radio spectral index between 408 and 1460 
MHz given by Vigotti et al. (1989). For the flux densities at 2400 A the fc-correction adopted was C = — log(l + 2) 



detail the relations between auv, z, P408 and L2400 reported 
in sections 6.4 and 6.5. The lack of an auv — P408 trend is 
still acceptable, in spite of an auv — £2400 correlation, since 
the latter is weak. The trends P408 — z and L2400 — z are 
likely related through the radio-optical correlation. 

If the correlation L2400 — z was predominantly intrinsic, 
rather than due to selection effects, the three correlations 
1/2400 — z, I/2400 — auv and auv — 2 could be explained with 
AD models with constant black hole mass and rh increasing 
with redshift (see Section 6.5). The correlation L2400 — P408 
would then imply that the P408 — z trend is, at least in part, 
intrinsic. If the trend P408 — z was due only to a selection 
effect, preventing the detection of low power sources at high 
z, the trend L2400 — z could be also the result of this effect, 
rather than cosmic evolution. In this case one of the cor- 
relations auv -2 or auv — £2400 could be induced by the 
other, in combination with the L2400 — z trend. According 
to the models by Wandel & Petrosian (1988), a correlation 
i/2400 — auv in the observed sense would arise naturally if, 
keeping a constant black hole mass, the parameter rh varies. 
We cannot exclude however the reverse interpretation, i.e. 
that the intrinsic trend is for auv — z and auv is not di- 
rectly dependent on blue/UV luminosity or radio power. 
O'Brien et al. (1988) found from a joint regression analy- 
sis for the three parameters that the dominant correlation 
for their data was between auv and z. 



8 CONCLUSIONS 

In this work we present optical photometry of a sample of 
radio quasars in the redshift range z = 0.4 — 2.8, around 80 
per cent complete, aimed at studying their spectral energy 
distribution in the blue/UV range. The U — R, B — R and 
U — V colours do not vary substantially with redshift, and 
the average values are (U — R) = —0.08 with rms 0.45, 



(B - R) = 0.64 with rms 0.27 and (U - V) = -0.38 with 
rms 0.42. Two quasars at 2=0.50 and 2=1.12 stand out for 
being particularly red, with U — R > 1. 

Power-law fits to the SEDs of the quasars with available 
photometry in the four bands yield spectral indices ranging 
from ~ 0.4 to ~ —1.7. The distribution of slopes is asym- 
metric, with a tail to steep spectral indices. Excluding the 
sources in the tail, the distribution is well modelled as a 
gaussian with mean and dispersion of -0.21 and 0.34 re- 
spectively. 

Composite SEDs nomalized at various wavelengths were 
constructed from the sample (exluding the two red quasars) , 
and the overall shape of the composites was found to be very 
similar for the different normalizations. The only emission 
feature revealed from the composites was the CivA1549 line. 
This result is in agreement with the expectations from the 
EW measurements of broad emission lines and Fen bumps of 
steep-spectrum radio quasars, which predict, for the band- 
widths and redshifts in our work, the largest contribution for 
the Civ line. For other emission features, like Mgn, Cm], and 
the Fen bumps at 2250-2650 and 3100-3800 A, the expected 
average contributions to the broad band fluxes are below ten 
per cent. The composite SEDs show a clear steepening to- 
wards high frequencies at around 3000 A which cannot be 
explained by line contamination and most likely reflect a 
trend of the continuum. Parameterizing the SEDs as power 
laws we obtained an average (a b iu e ) = 0.11 ±0.16 for the 
range 4500-3000 A and {auv} = -0.66 ±0.15 for 2600-1700 

A. 

Separating the quasar sample in two redshift bins, with 
the cut at 2 = 1.2, a better agreement was found between the 
values of abiue and auv obtained for the different normaliza- 
tions, and a hardening of auv with redshift emerged, with 
(a UV ) = -0.87 ±0.20 for 2 < 1.2 and (a UV ) = -0.48 ±0.12 
for 2 > 1.2. The average spectral index abi uc for the low 
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redshift quasars is (abiue} — 0.21 ± 0.16, therefore these 
quasars show a steepening from the blue to the UV range 
(abiue) - (auv) = 1.08 ± 0.26. The composite SEDs for the 
high redshift quasars do not cover the region above 3000 A, 
therefore the presence of a similar break for these quasars 
could not be analysed from the present data. 

Correlations were also found between luminosity at 2400 
A and redshift and between radio power and redshift. Sepa- 
rating the quasar sample in two bins of low and high lumi- 
nosity at 2400 A a trend was found between auv and L2400 
(89 per cent significant). The quasar sample shows also an 
intrinsic correlation between L2400 and P408 (r s = 0.62 and 
P > 99.9999), similar to that recently reported by Serjeant 
ct al. (1998) for a sample of steep-spectrum quasars. 

The observed trends L2400 — «uv, £2400 — z and auv — z 
appear to be consistent with the predictions from AD mod- 
els for the case of constant black hole mass and accretion 
rate increasing with z (Wandel & Petrosian 1988), with 
Mbh — 1O 8,2 ~ 8,5 M0, and accretion rates (in units of the 
Eddington value) ranging from rh ~ 1 for z = 0.6 — 1.2 
to rh ~ 10 for z = 1.2 — 1.9. An alternative interpretation 
is that the P408 — z and L2400 — z trends arise predomi- 
nantly from a selection effect, due to the radio flux limits of 
the sample. In this case one of the correlations, auv — 2 or 
auv — 1/2400 , could be induced by the other, in combination 
with the L2400 — z trend. The observed auv — £2400 correla- 
tion is consistent with the predictions from the AD models 
by Wandel & Petrosian assuming a constant black hole mass 
and a range of accretion rates, although in this case the ac- 
cretion rate does not need to be correlated with z. We cannot 
exclude however that the intrinsic correlation is a\jv~z, with 
auv not being physically related to the blue/UV luminosity. 
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